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We present updated results for the radiative improvement of the a ■ B term and the spin-dependent 
four-fermion terms in the lattice NRQCD action, and first results for the radiative corrections to 
the NRQCD Darwin term and spin-independent four-fermion terms. The spin-dependent terms 
have significant impact on getting the correct hyperfine splitting for both bottomonium and heavy- 
light mesons, while the spin-independent terms suffer from a conspiracy between lattice artifacts 
and severe IR divergences that complicates their evaluation. 
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1. Motivation 

Non-Relativistic QCD (NRQCD) [1] has been used successfully to describe both quarkonia 
and B q mesons. However, until recently, only the tree-level NRQCD action has been used. What is 
now known is that the radiative improvement of the coefficient of the a • B operator in the NRQCD 
action and the inclusion of spin-dependent four-fermion operators into the action have a significant 
effect on the value of the bottomonium hyperfine splitting [2]. Here we present an update on 
the improvement of the a • B and four-fermion operators, as well as first results for the radiative 
improvement of the Darwin term in the NRQCD action. 



2. Matching NRQCD to QCD in Background Field Gauge 

The NRQCD action used by the HPQCD collaboration is 

S = £ ¥ f (x, t) [y(x, t) - K(z) y(x, z)\ 

.V.T 

with the kernel 
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and n > 3/(Mo«) is a stability parameter used to avoid the well-known instability of the Euclidean 
Schrodinger (i.e. diffusion) equation. The normalization of the operators is chosen such that the 
tree-level matching ofNRQCD to QCD (which can be performed e.g. by a Foldy-Wouthuysen-Tani 
transformation) gives c,- = l + ff(a s ). 

In order to match NRQCD to QCD at the one-loop level, we have to determine the radiative 
corrections to c, by demanding that some suitably chosen set of renormalized S -matrix elements 
agree to one-loop accuracy when calculated in QCD and NRQCD. For the quark bilinear terms 
in 8H, we can use quark scattering off a background field for this purpose. Hence we are led to 
demand that QCD and NRQCD give the same effective potential after non-relativistic reduction, 
i.e. that the following diagram commutes: 



qqj2 NR reduction tree-level 
using c, 



1PI 



NRQCD 

1PI (2.1) 



NR reduction 
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2.1 The Background Field Method 

The effective potential in the presence of a classical background field <I> is defined by 



where the path integral over the quantum fluctuations only makes sense in perturbation theory 
restricted to IPI diagrams. 

In a gauge theory, we can decompose the gauge potential into a background and a quantum part 
as Aji = +gqfi. BRST invariance of the classical action S guarantees that all D < 4 operators 
appealing in the effective action T are gauge covariant, implying the renormalizability of the theory. 
The D > 4 operators that will appear in T are, however, not necessarily gauge covariant. 

As an effective theory, NRQCD contains D > 4 operators in the classical action S, where gauge 
covariance can be imposed at tree level. At the loop level, gauge covariance must also be retained to 
avoid serious complications (even though the appearance of gauge-noncovariant operators would 
not render the theory invalid in and of itself, keeping track of their gauge dependence would be 
highly cumbersome). As has been shown long ago, [3], this can be achieved by using background 
field gauge (BFG). 

Background field gauge is defined by the gauge fixing function 



and hence not only the propagator, but also the three- and four-gluon vertices of the form qqB and 
qqBB are gauge-parameter dependent. 

An important feature of BFG which we will exploit in the following is that it leads to QED-like 
Ward identities and finite counterterms, so that we can compute all diagrams numerically, and do 
not need to calculate the gauge field renormalization. This is practical in particular for checking 
the gauge-parameter independence of our results (note that c\ and C2 should be gauge-parameter 
independent, since they are directly related to physical mass splittings). 

To employ the background field method on the lattice, the gauge link is decomposed into the 
ordered product 



leading to a dependence of the Feynman rules on the number of background and quantum fields 
(qqq, Bqq, BBq, etc.), as well as to the appearance of different terms for different orderings (Bqq, 
qBq, qqB, etc.) contributing to the same vertex. Background field gauge (BFG) is defined by the 
gauge fixing function 



which on the lattice affects all vertices with exactly two quantum gluons. 

The automated derivation of Feynman rules with background fields, as well as the use of BFG, 
has been implemented in the HiPPY/HPsrc packages for automated lattice perturbation theory [4]. 

Lattice gauge theories in BFG are renormalizable [5]. Checking that the gauge dependences of 
our results match for individual terms, and that their sum is non-trivially gauge independent, gives 
us confidence in the correctness of the results. 
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Figure 1: The one -particle irreducible diagrams entering the matching of the a B and Darwin terms in 
NRQCD at the one-loop level. On the QCD side, only diagrams (a) and (b) need to be computed, while 
in NRQCD all diagrams contribute; note that diagrams (c)-(f) contain not only lattice artifacts, but also 
physical contributions from higher orders in \/m. 

2.2 Matching the spin and Darwin terms 

The effective action for continuum QCD contains terms of the form 

rpp,y,A] =w 1 ( g 2 )^y+y 2M ^ v ^ + --- 

which after renormalization and non-relativistic reduction give 

(l+F 2 (0)F 1 (0)- 1 )Y^YR 

n v ZMr 

for the a ■ B term in T NR , and 

(1 - 8MlF/(0) + 2F 2 (0) )v4 ("f|r) Vfe 

=b D 

for the Darwin term. A straightforward continuum calculation gives 

l( i) / 3 , n 13 \ ,m ( M 1 1M 1 50, ix \ 

On the other hand, the effective action for NRQCD contains the terms 

m , a] = C4 zS v v - <*% V ^ v + ... 
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which after renormalization give 

T^hir i/f 1 " 41— ^ - 7 NR- 7 NRyNR 1 ,.t ' a ' (g „, _, y NR yNR / yNR \ 2 . , ,t ggj^O , 

1 LVr> Vr> a J - c 4Ar z 2 z m % — ^j~ R — y R -c 2 z, D z 2 (z m ) Yr gM 2 Vr + --- 
Equating the corresponding terms in T NR gives the matching conditions 



„ V NR V NR V NR i i . 

c 4 Z a Z 2 Z m = l+b a 

„ vNR v NR,"7NR\2 i , • 

c 2 Z D Z 2 (Z m ) = \+b D 

M 



which yield c\ — 1 at tree-level, and at one-loop order (with Z = 1 + 8Z) 

c? = b$ - 5Z^' (1) - 5Z 2 NR ' (1) - 25Z ; r (1) 



In lattice NRQCD, we also need to take into account contributions from mean-field improvement 
U i— > U /uq in SZ^yj^ besides the diagrammatic contributions. 

2.3 Four-Fermion Operators 

Beyond tree level, the NRQCD action also contains four-fermion terms 

which can be rearranged by a Fierz transformations into 

2 2 

= a l jf f2 (^x)(Y t Y)+a^(x' l t t a x)(Y l taY) 

2 2 

+b l ^(x t °*X)(Y t °Y)+b 8 -^(x' f o*t t a x){y Jf ot cl y) 

where the coefficients d{ are linear combinations of the coefficients a,, bu with the latter directly 
computable from the box diagrams of fig. 2. 

In addition to the box diagrams, there are additional contributions to d\ from QQ annihilation, 
which is possible in QCD, but not in NRQCD. For our purposes, the relevant contribution is 

In 

d r = ~(2-2lo g 2). 

3. Divergences and artifacts 

The calculated quantities in both QCD and NRQCD are IR divergent. We can use a gluon mass 
H to regulate these divergences, since in the absence of ghost loops, all gluon lines are attached to 
a conserved current removing the unphysical longitudinal mode contributions. The NRQCD con- 
tribution then contains an IR logarithm log(/xa), which combines with the IR logarithm log(;U/M) 
from the QCD side to give the expected logarithmic log(Ma)-dependence for the matching coeffi- 
cient. 
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Figure 2: The one-particle irreducible diagrams entering the matching of the four-fermion interactions in 
NRQCD. In QCD, only the box and cross-box diagrams contribute, while in NRQCD all diagrams must be 
computed. 



Figure 3: A graphical depiction of the lattice artifact contribution from the insertion of the gluon self-energy 
on a single Coulomb exchange line. 

On the other hand, any power divergences must match between QCD and NRQCD. The spin- 
independent part contains power IR divergences up to M 3 //i 3 , whereas the spin-dependent part is 
much more mildly divergent, with IR divergences only up to Mj\i. We treat the power divergences 
by subtracting the known divergences analytically from the NRQCD integrands before performing 
the loop integrations. 

This is enough to make the spin-dependent part well-defined. In the spin-independent part, 
however, the presence of k 2 a 2 lattice artifacts can move the leading IR divergences to lower order, 
leaving an apparently meaningless (Ma) 2 log(jua) overall divergence in the sum of the box dia- 
grams. The origin of this artifact divergence can be seen to be from the octet Coulomb exchange; it 
cancels against the artifact divergence arising from inserting the self-energy a 1 correction on a sin- 
gle Coulomb exchange line, which also contributes to the lattice artifacts of the spin-independent 
four-fermion term. Since the spin-independent part is numerically very difficult to evaluate due to 
the presence of these artifact divergences, the results for the spin-independent part are still under 
review. 

4. Physical Impact 

The bottomonium hyperfine splitting (HFS) receives a contribution ~ c\ from single-gluon 
exchange, which is corrected by the spin-dependent four-fermion interaction ~ [d\ — d%). Empir- 
ically, the latter contribution is found to reduce the lattice-spacing dependence; a larger value of 
the four-fermion coefficients also partially compensates for the smaller HFS found when including 
spin-dependent ^(v 6 ) terms [6]. The physical impact of the radiation improvement is significant: 
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the measured IS HFS ~ 60 MeV is corrected to 70 MeV, in agreement with experiment [7]. In 
addition, the HPQCD prediction [8] of 35(3) (1) MeV for the 2S HFS including the radiative cor- 
rections agrees with the newest Belle results [9]. 

In the heavy-light case, the HFS ~ C4 receives no four-fermion contribution. HPQCD finds 
good agreement with experiment for the Bj, B s HFS when including the radiative corrections to C4, 
and is able to make a prediction for the B c HFS [10]. 

The impact of the Darwin term is much less noticeable: the bottomonium S-wave energy shift 
~ C2 constitutes a very small effect except on very coarse lattices. 

5. Conclusions 

We have extended the radiative improvement of a • B term and spin-dependent four-fermion 
terms in the NRQCD action to the n = 4 v 4 and v 6 actions, and have computed the radiative im- 
provement of the NRQCD Darwin term for the n = 4 v 4 action. We are currently computing the 
Darwin term of the v 6 action, as well as the spin-independent four-fermion terms. 

The spin-dependent terms have a significant physical impact: in fact, the agreement of the 
theory with experiment depends on radiative improvement. The spin-independent terms are more 
subtle to compute, but have much smaller effects on heavy-quark spectra. 
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